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FOREWORD 


eepite tresendous technolojical discoveries and developments, 
the gas and ail inilustry of today is still plagued with unresolved and 
unexplained problems an) phenenwsna, ospecially in the field of reservoir 
engineering. me of these is the sacrescupic flow of fluids in reservoir 
reek formetiona where «as and oil are found and sreduced, This thesis is 
a stady ef but a email facet of the creblea and has been conducted with 
the idea that any work dome in this ares will heip contribute to an in- 
creased knowledge of the subject. 

the auther wishes to seknosleize the simere end helpfal eugces- 
tions, advice, and encouragenent of Profeseor “olbrook ¢. Rotest, ‘ead of 
the Patroleus Engineering Seperteent, University of Pittsburgh, in the 
projection, progression, and suecossful eanpletion of this stuty. 
Aeknowledgment ia also made to those authors, investicatora, sid researchers 
whose works were utilized for background end basis of the study without 
which such en undertaking would have been iepossible, 
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4 . LRTZING TT ON 


The Patroleunm Industry finia and obtains ite ras waterials, cas 
and oll, far beneath the sarface of the earth in areas of rock formation 
ealled reservoirs, These reservoire are not, as is usually presosed by 
the uninformed, large voide filled with gee and oil which, when tapped, 
gay be eaaily produced. Om the contrary, they are formations of porous 
and cormeable rock with small, wacroscepic, open spaces between the 
grains of the satrix saterial in whieh the gas and oil is stared; a 
fagtor which makes the croduction of such stored gas ami o11 difficult 
amd, in some ceases, ianossible., itis for this reesen thet the coapiete 
study of all comiltione an4 probabilities of fluid flew in porous sedis 
becomes an incertent facet in the preduction of gas and 011, 

The purpese of thie study was to investizate and attennt to 
evaluate the inter~-relationships of some of the major physical proverties 
of porous reservoir rock which affect the flow of Muids through it. 
Although these inter-relationships have been expressed by several di f- 
ferent formulations, mentioned later, 612 designed te elucidete fluid 
flow, it was thought that a mere useful and informative modi fication 
could be developad ami applied to the problem of correlating permeability, 
poros!ty end grain aise of porous wadia to fluid flow through such media, 
Sueh a atady, designed to reveal the wsechenics of secrescopic fluid flew 
in porous media, shoalé be usefal in evaluating some of the sany factors 
whieh mast ultimately be known if a satisfactory solution of the croeblea 
is ever to be obtained. Any adiitional knowledse of these factors should 
contribute to a more conplete understanding of the conplexities of cas 
and oil sroduction from the reservoir, 
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Yo the petroleum engineer, especially ome interested in reservoir 
engimering, @ very important. physical property of reservoir rock is 
permeability, which can be properiy defined as a massure of the fluid 
transeiiting eapacity of & porous material, or ag the abllity of a fluid 
to flow witnin the iaterconnected pere spaces of a poroun material. A 
aecom! physical procserty, also equally important, is poresity. In ofl and 
gee reservoirs, porosity represents the pereentage of total reck volume 
whieh is avallable for occupancy by either liquids or gasses or both. In 
connection with the term porosity, care should be taken to distinguish 
veterwen absolute and effective porosity, Absclute poresity is the per- 
centage of void espace with respect to totel volume. Effective porosity, 
on the other band, is ‘he pereentage of interconnected void apece in the 
total voluge, A third physical property of reservoir rock is actual 
grein size eni grein size distribution of the matrix saterial which forns 
the reservoir rock. 

lt has been demonstrated by experiment and empirical fornmlas 
that permeabllity, effective porusity, and grain sive are interrelated, 
one depending upon the other te varying degrees. In porous material in 
which the pores are inter-connected, there exiate no definite relation 
between permeability ami effective poroceity. ‘These, however, are ilrectly 
influenced vy grain size, in that grain size or particle seize aw! ite 
attendant, fectors of uniforsity of grain sise, shape of the crain, sanner 
of grain packing, and the amount of cementing material between the grains 
(degree of lithifleation) deternine the size and shape of the pore open- 
ings, the extent of intercorrections of the pore spaces, and the ratio 
of their volume to the total volume of the rock. 
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have indleated that actual grein else 
and effective porosity are excellent criterla for judging permeabilities 
of various porous media, in mach of the literature, variations of 
peraeability and porosity within iniividual reservolre have been found to 
be inter-related to a degree and have been shown by Fettkels? in the case 
of tra Bradford Third sand of the Brecferd field of Permmsylvania te have 
these relationshipa 

K = 6eu9 x 107% g )5e85 
and K a Sad x 197204 g )°*? 
where K is the permeability and J is the perosity, The firet formule vas 
deaterained from earlier exceriwents and the second from later experiments, 
However, thase milationships are average and do not apply in some cases or 
taroughout the complete possible range of perseabi lity and porosity, ere 
again grein sise, with ite attendant factora, directly affects these 
relationships. 

In waking deteruinations of the values of orrmeatility, effective 
porosity and orein size diaseter, the first two are now very easily and 
effectively carried out by meang of slaplified, standardized metheds which 
give fairly acourate results, in the case of average grein size dieneters, 
thie detercination becomes complicated by the particle iteelf, in that the 
particle is ugually irregular in eise smd sbepe, and no know sothois are 
availeble for defining such a particle in geometric torma. Thie problem 
becowes core comlicated ani forsidable in the case of an aygregate of {r- 
regular particles, but has been successfully solved by the use of statiat! - 
cal methodea for evaluating surface area and volume of the aggrecates, 


1 References sre listed in the Mbiiography. 
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Such sothois are bewed om the aeasurement of the so-called avaraze crain 
ameter by rwiating the irregular particle to an equivalent ceomatrica 
figure such as a circle or a sphere from which a diameter can be calealsted. 
These saethods deteraine surfece area ani volume froa the number end 
weight of tne particles, using o@ a besia the theory of ciaensiens which 
requires that tae surface area be proportional to the square, and the 
volume be prosertional to the eube of the dimenelons used. 

For the deteraination of particle sian to have physicel signi fi- 
cence, that is, to relate particle size to ‘efinite physical properties 
sugh es surface or volume, the aritimetic or ceometric sean and the netian 
dlLaseter ceterminations are supersede: by statliatical diunmsterm. Zeveral 
formulas have been developed and used”? to designate these statistical 
diaseters in terw of various physical properties such as acan dieneter 
(length), meen volume diameter, moan surface ‘lameter, moan volune-surface 
diameter, end weight sean diametes, The average grain dlemeters ueed in 
this study are mean volume diameters, based on the average volume and de- 
fined to be that diameter whose volume (divided into the total volume would 
clive the number of particles. If dis the average particle diexeter and 
¥ ia the totel volume, & relationship existe es 

v aS nd? 
where nis the nusher of particles of a given size, If 4,? represents the 
volume of an average particle, and & nis the total number of particles 
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Tle ARSTORICAL RACK ROUND 


Investigation of the flow of fluiie through porous asterial at low 
pressures wos imtiated by -arey? in 1556, ‘These atuties consisted of a 
series of experiments on the flow of water through filter bede frea which 
varey obtained an explrical foraulea showing the rate of flow to be very 
nearly proportional toe the pressure drop per unit lensth of porous medium. 
In lte eloaentary fora, this formla, 

va ~(K/1) OPAL) 
iniicates thet the fluid permeabliitty, K, of a porous material proportion 
ately relates the velocity, v, of flow of a fluid of » viseoslty to the 
preacure differential, 4S), causing thie flow over a length of OL. Since 
QA? is asasared between the outflow and inflow ends, it is negative as in- 
dicated in the formule above. 

The propertionality constent K is considered to be a specific 
property of the porous material, eapirically iniepencent of the dimensions 
of the auterial, the pressure differential exerted on the fluid flowing, 
and the viscosity of the fluid. For this reagen, K may be expressed in 
terme of other measurable physical properties of the porous material, such 
as porosity and grain aise. Cince the formulstion of lerey's equation and 
fron it, his law, uany investigaters and researchers heve attempted to show 
relationships of porosity, permeability, and crain size, 

in 1800, Seeiheist! first introduced grain size into the relation- 
ship by showing that the rate of flow was proportionsl to the squere of 
the average crain diemeter. ‘This would indicate a decrease in permeability 
aa sand graine decrease in size. 
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the next investirstor to utilize yrein size in flow relationship 
formless was teson” who, in 1592 snd 1599, proposed the formla based on 
studies of sand Miltere, 

q « &,, ag) (n/t) 

where K. ie @ permeability constant, h/L iw the pressure drop through the 
filter per unit=iepth, and ay, Le the "affective" diameter of sand particles 
havin; aniforulty coefficienta of less then 5. "iffeetive" crain size is 
defined es the epentng which wili just sess ten per oent of the purticles 
(ty weight) end uni foraity coefficient 1e sefined as the ratio of the sise 
opening which will pass 60 per cent of a sanple bein: screened, to the 
sise which will juet pass 16 per cent. Ehen the uniforsity coefficient ‘is 
low, ea 5 oF less mentioned above, the particles are sore or leas unifons 
in size, but as the uniforsity coefficient inercases to higher values, the 
particle sixes becom widely distributed. 

in u@ing fasen's equation, if q is expreased in cubic feet of 
water per day per squere foot of filter srea, 4, in willimeters, the loss 
of head, h, in foet of water, ani Lin feet, thea the veluce of K,, range 
from 1300 te sebove i000 descending on the cleaniiness of the semi, with 
the usual limits of K.. fer ordinary sande ranging between 2300 ami 3300. 
Hagen's equation has been widely accepted and utilised oy sanitary ongincers 
and is of considerable historical interest since 1t represents the first 
attespt to recagnise the importance of partiele size and a aethed for ite 
satisfactory representetion. Savia and Filsey? retested Nasen'a Zormla 
with the following resalte: 


as 2300 4,2 (ELA) nj) 


, 6 [2¢,39) ,, ». 
- @ = 780(4y,)*(#/u0) ox 
where T is the tesperature of the water in degrees Fahrenheit, dy, Le the 
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screen sise in alliiasters through which just 34 per cent of the cartleles 
{ny weicht) pass, and g is the porosity inixedueed inte the second formals, 
Sidenter,”’ in 1899, published the aecounts of hie theoratical 
analysis of fluid flow through an ideal homegeneous pecking conpoeed of 
apheres of uniform si.se in which he was tha first to introduce the effect 
of packing a3 @ factor influencing pereeability, which <ave rise to ex 
pressing the average pore grea in terse of the diameter of the sphore,. 
Slichter cerived a motified Poleeullie relationship in his equation 
qs l0.2hd*h/K ul 
where q ig expressed in ac of water per secom!, A ie the crosa<sec tional 
area of the pecking in square centimeters, h ia the presoure differential 
in centimeters of water, d is the diameter of the spheres in centimeters, 
K, is a packing sonstant cepenient on porosity, 4 is the visessity of the 
fluid flowing in poises, anc L is the thickness alon;s the direction of 
flow in centimeters. in teness of curcy's permeability, Slichter's formala 
would read 
K = 10.24°/K, 
For various packings, tabulations have been nade of the different values 
of the pecking constant k, shown ag a function of poresi ty fer perous 
eaterial packed with spheres of «qual claweter. 
Zavis ari uileey” Seveloped a formla 
I/i, = 0005(6/00)?*? 
for use in the determination of the values of Ky where # is the peresity. 
ting’ recognise’ the importance of Slichter's formule by using it te ecal- 
culate the averace ‘Lemeter of irregular particles, It is importent to 
note here that the average dlamater of irregular particles determined in 
this menner 1s the diemeter of the perticles such that if a1) of them were 
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of this diameter, the vaeking woul! heve the seme permeability es it 
actually hes at g civen temperature and poreslty,. 

For the ultimate in permeability, porosity, and crein sise rele- 
tionships, it remained for keseny,° in 1927, te develep and publish his 
ispertant empirical formule 

v a (Pe OP/a8*L) (u/ty) 
where C is a form facter dependent on the shape of the pores, § is porosity, 
g ia the aeceleratien constent, Oi’ ie the pressure differential in gram 
por square centimeter, u ia the viscosity of the fluid Mowing in poises, 
5S is the specifie gurfece area in square centiseters per cuble cartineter, 
Lie the thickness of the porous nediuse in centiaeters, L, is the longer 
path of flew in centiesters through the thickness L, and L/L, ie a redue- 
tion factor for the relationship 4P/i. Since § @ 5,(1 - 9) the above 
eqeation cen be rewritten in temas of S., whieh is the specific surface 
area per unit volume of setrix material, ea, ae 

Ya (oPg ar/us* 2 ~ £)*L)(L/t,) 

Unaware that Eoseny had develeped the afcrenentionsd relationahig, 
fair and Hateh,”? in 1933, developed an almcatidenticel formla by folles- 
ing @ line of reesonin: similar to that ourgsed by Koreny. This equation 
is | 

¥ 2 PaP/Sall « fy 
It ean be ween that the fastore (C*L/i_) in the Koseny equation have 
been replaced by the constant 5, believed by Fair and fateh to apply to 
normal, rather compact, unconsolidated porous asterial, 

CarmenOrtl in 1937, presented an improved derivation of the 
Loweny ani Fair and Hatch equations, Using ea @ basic the gereral equs- 
tion for flow throagh pipes 


08 GS WILE es ce OP ote) Cem ee we tere tk * 
Ole, oe eres Bea) gM ote Clete 
ple Me le, ee ee 
Sey ee of@ec at *, gree 5? ree 8. ee li 6 
' ee ee hee 


7“ 7™ ee ee Tee 





NE AEN yee a Se eg Oe ae pO were anti A no btm 
Se 2 CET re oo ot Aker materiince: Ab 4: 5 
sRonkny 0) Atm) Sd? oe Ve RhOmEEEY ue AS + yrerm yet moe IG 

ST moge se ait 





mee Oke Saad ea er les eabietin te ta ‘a 
a 
* Aneta tage 


ores mimes yates att cs “WO emt ae ee ee we et) 
4 cil St Ay ae A ow Semrdins 44 rede ob @ bent: amet 
aly ‘aris Wall a iaeee sirciiiaeas eiian samaiin 
eee = 
aw Sas ass ee ANS GET ol tn eae 
fh ett ea ek ant 





v= gn 4 P/E, 
‘whe De : is the average velocity, ¢ is the acceloration constant, ais the 
hydrwalic radios, defined as the volume of fluid in the pipe per surface 
area in contact with the fluid, A? ie the pressure differential, k, is a 
shape factor ami is the reciprecal of Koeveny's form factor ©, u is the 
viseosity of the fluid flowing, end L, le the length of flow path, Carmen 
set about to substitute or movlty the equation facto7s so that the equa- 
tion would represent flow of fluid through porous media. For v he sube- 
stituted v/f, later to sodify it by the relationship L,/l fer flew over 
the flew path, l,, creater then the i dimension of the medium, Then he 
substituted for a the relationship of $/5,(1 - ¢), and applying the redue~ 
ing factor L/L, to the expression 4 ?/L, be arrived at the following 
rearranjenent : 

vs 00/6. 2a = #)*LNu/y) 

in a series of experiments using @ very wide range of particle 

shapes in unconsoliiated media, Camsan found that the peraveter (CeL/i,)~> 
of Koseny'a equation and the parameter ka(L/l.)* of his equation both 
averaged sbout 5. In the moiified oseny~Carman equation these parameters 
are expreseed as ky with en assigned value of 5, shown by 

ve A/a kt ~ #)*L 

#yllie, in his paper on the historical development of the 

Kosery—Carman equation, discusses the sere recent detercinations by ine 
vestizators of values greater than 5 for the Kozeny~Carsan constent k 
when the formla hes been used in connection with measurementa of flow 
through consolidated maiia, Since the k, shape factor of the Keseny~carsuan 
equation has beon held by scost investi,zators to be constant, the hi cher 
values of k have been attributed to the (14/1) factor which fose end 
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Brace’? have called the tortuosity factor, T. Several determinations 
were made hy ose end Bruce for the value of & for various consolidated 
aodia of low permeability where values were obtained as high as 100. 
Others concur with Soese and Bruce, ali iniicatin; that the constant for 
consolidated matia varies to higher values and that these hicher values 
are attributable to the tertucsity fanction. 

Thers has been some disazreement with the concept of increasing 
values of the Koseny-Ccarsan constant in connection with censolidated media, 
le Ee Fose™® pelieves that errors arising froa the use of 5 as u value for 
the constant will be far loss than those arieing from uncertainties in 
other variables, especially the poresity function, expreseed as P/(1 - g)* 
in the Koveny~Carman equation. fie hes shown, by plotting relative 
resistance versua per cent porosity, that at porosities lower than 4&O per 
cent there are increeaingly <cifferent values of relative resistance to 
flow for various equations which all contein a porosity function of som 
sort. oabLevalie®” hes shown @ similar plot of relative permeabliity versus 
per cent porosity for several different inveatigators' fornulss, all con- 
taining « porosity function. This algo inileated that st porosities lower 
than LO per cent the porosity function does not remaln constant. It has 
net been possible to test the porosity function 9/(1 - ¢)* in consolidated 
nedia because (t ls not reedily poeasible te vary perosity and keep the 
specifie surface area constant. 

A vecent moi fication hes been made in the Koteny equation by 
fylide and Spangler! who have cosbined the Koveny equation with properties 
of the capillary pressure deasturetion curve to obtain the formula 

Kw (v7/2.5776) 5 a5 /?,° 


where K is the permeability coefficient, y le the interfacial tension, F 
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is the formtion feetor, ¢ is the porosity fraction, S, ia the fractional 
wetting phase saturation, and pe i3 the eaplllery pres@ure. This equation 
introduces new physical factors wiich are more easuily neseureble in 
leboratory jeterminations and consequently should further the continuing 
eearch for nore accurate and useful reletienship formulations for fluid 
flow threugh porgus madia, 
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ili, PROCERORE 


For the purpose of carrying out thia study, data cof experimental 
and calculated values for persenbility, together with secompanying values 
for porosity and average grain ise, wore selected from Fettkets'?* 
atucdies of the Bradford 911 Field, Bradford, Penneylvenia. These dete 
had been derived fros test deterninations mee on representative core 
sanplesa obtained from various sections of the fleld at random destha in 
the Bradford Third eand. The Krecferd Third sand is on excellent rezre- 
eentative of a consolidated porous wedium whose seneral characteristics 
are low permeability and porcelty and quite fine créin eise with irregular 
shape and aise distribution and a moderate to high ancunt of 11 tht Mestion,. 
Sines practically all cil ani gan reservoirs are lecated in consolidated 
roek formations, a representative sancle of tiia type of porous setius 
was chosen for application to this study. 

4¢ @ first approach te the problem, an investigation wes made of 
permesbility, porosity and grein siae relationships when applied to an 
cuuiies wetithet wy the eather tem om teva ey Teber ent be” 
from the Poiseullle equation of flew threugh pipea. This equation wee ex 
pressed wy Traxler and Baus as: 
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Ne = (32qul./ag oP) (2) 


where D, ia the effective pore diameter in centimeters, q is the rate of 
flow in cubie centimeters per second, j. is the viscosity ef the fluid 
flowing in peises, L ia the length ef flow path in centimeters, A is the 
crose~sectional area in square centimeters, ff ie the porosity, and AP is 
the pressure differential in dynes per square centimster, It can be 
readily seen that the function q:l/AO?P of this equation will be equivalent 
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to the perseability coefficient K, in darcies, when the units are abange 
to those of the Darey equation, Thus, K a (G100sL/49.8692 ¥ 1077 AP) or 
95692 X 10°? EC 2 @L/AAP, Sudetituting in equation (1) above gives: 


Dg = 5629 x OKs)? (2) 


where K is the permeability coefficient in dareies. 

Table I, Appendix I, has been presented to show calculations wiich 
resuited from the aprlicsation of the permeability and porcsity dete selected 
for this study to equation (2) abeve, solving for Dy, the effective 
averaze pore dianster. Hutting”” hea stated that the effective averege 
pore disaster is roughly one-fifth of the average grein dlanster, Pigure i, 
Appendix II, shows « plot unde fron the data in Table I of d, the average 
grain dlemeter, versus Oye the effeative aeveregs pore diameter, resulting 
in a straight line whose slope, which represents the ratio of d/),, 1s 
approximately 25. The value of this retio is considerably crester then 
the one of 5, (1/Sd 2 ,), proposed by Rutting, It should be explained 
here that enly the rweulte of the caloulationa using Fettke's 1934 date 
were used in plotting Figure 1. The results of the caleulations using 
Fettke's 1933 data praiueed no conclusive plot. 

AS & Second approseh to the problen, an investigation was made of 
permeability, poresity and crain sise relationships when applied te a 
formulation derived by equating the Koseny-Carman and the Carey equations, 
This second study was made on the dagis of the specific surface area, 

Sq, of porous medium per unit volume of matrix material, Sy taking the 
basie Koveny-Carsen equation solved for 5,” one obteins: 


2 ger. ar (3) 
S(l = #)* qu 
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it. can be sean thet in equation (3) om the preceding page the function 
(AA P/qab)7> Ls equivalent to the Ko of Carey's ecactien when the write of 
} and O P era chanysed from polses to centipolees an! from grams per 
square centineter to atmospheres respectively. Thus, K = gl0GuL/ 


=< 


(92673 x 104?) or 9.678 x 10 £ = (AN P/qul) >. Substituting in 


equation (3) zives: 
2 


3 
aa” (iy) 
o = Fa - i.e x 1 


fince the acceleration constant, g, is 930 centiaaters per second equered, 
equation (4) above can now te resolved as 


Sq = USCls ar (o/iy’ (s) 
woieh is the equation used in this seconi ef reach to the oreblem anter 


atudy. 

Table Ii, Acpendix I, hea been prepared to ahow calculations re~- 
sulting from the application to equation (5) above of the permeability 
avi porosity data selected for this study. 13 tee bom Game & Suen 
that S, = Z,/2,4 for irregular particles where %,/Zy 19 the ratio of the 
surface to volume factors end d is the average crain iieseter. In this 
ratio, I, 2 G/d,", where S is the ourfoce eres ani d, is the statistical 
particle diameter based on surface arra, ani 2. zs v/ty?, where Y is the 
volume and a, in the statistical particle diameter based on volume. ‘The 
ratio Z,/é, has been iniieated ty Fair and fateh” to be a useful measare 
of particle shape am! has been <iven @ value of 6.0 for spheres, 6.1 for 
rounded particles, 6.) for wern particles, 7.0 fer sharp particles, and 
7? for angular particles. Ficure 2, Appendix Ii, shows a plot made cf 
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the values in Table li of &, versus 1/d whieh roughly for: straight lines 
with & slope of 143 for line 1 wailing Fettke’s 193) data and a alope of $7 
for line 2 using Fettke's 1935 data. The slope norsally should represent 
the values of the ratio 2/2, which lie somewhere between 6,0 and 7.7. 

Aa @ third approsch to the problem, an investigation was made of 
pormeability, pemsosity and grein aise relationships when applied te a 
formation Jerlved by equating the Aoseny-Carman ani the “arcy equations, 
a8 previously mentioned, except that the formiation wes modified for the 
use of average grain size in lieu of specific surface area, Since 5, is 
proportional to 1/d, s,* can be replaced by v/ae in equation (1) above. 
There is also a new constant, k', introduced into the equation whieh Le 
equivalent to the function IW (Eg/ty)”S where S is the Koveny~CSarsan con= 
stant and 2,/Z, is the ratio of the surface to volume facters. The re- 


gulting equation is expressed as: 


3 
= 2 TT. (6) 
(1 = §)?9.673 x 20 


By oudstitutin; 900 sentimeters per second squared fer g in equation (4) 
end rearranging to solve for k', the equaticn can now be expressed aa: 





xt = x (a - g)* (7) 
d 1.013 x1 

Table Til, Appemiix i, hae been prepared from caleulations deri ved 
by applying the data used in this study to the fectors of equation (7); 
namely, K/é® ana 1,013 x 10°/(1 = ¢)", ‘the values of these factors were 
plotted ae x/d* versus 1.023 1 10°99 /(1 ~ g)*, shown by Figure 3, 
Appendix II, to be a atraight line with a slope of 0,.0001052 in the case 
of line 1 using Fettke’s 1934 data and a curve with os variable slope of 
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unknewn value in the cave of line 2 using Fettke's 1936 data, iine 2 
could not be expressed a8 a straight line by repletting on sent-loy or 
lerm~log graph paper, The slope of line 1 of Migure 3 is equal to the 
constant k* of equation (7) on the preceding page. 

he previously stated, k' is equivalent to the funetion 1/(2,/z,)°5. 
Uaing the value of the slope determined for line 1 of Figure 3 for k', the 
value of 2/2, as obtained from the relation Zg/ty = (1/k'5)* 48 4305. 
formally, values for z fe, lie between 6.0 and 7.7 a8 previously aentioned. 

In the two preceding investigations, the perosity function, 
g°/(1 = g)*, sithough 4t varies with the value of porosity, ia considered 
@ constant factor in that it retains the same fined relationship recgarid- 
lees of the source of data. For a fourth approach to the problem, an in- 
vestization was mede of the permesbllity, porosity and grain size relation- 
ahipewhen epplied to a formulation developed to expresa the porosity fune- 
tion aa @ varishble, dependent upen the source of data, with the perses~ 
bility and grain size belng constent factors regardless of the sources of 
date. ‘Thus, the porosity function would be s constant factor only for «a 
particular reservoir formation or even a particular section of the 
reservoir formation. 

To develop this formlaticn, the basic Kozeny-Caraan equation vas 
recast to include particle ileneter, giving: 

‘a * pp (8) 
a L (1 = 8) 

whose units have been previcusly defined. Solving for the porosity funetion 
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~ z a Le (9) 
ci | ged (re Be gie 
(1 = #) 


AOP z k' as 


Since the function qiL/AD >, when expressed in teres of the “eray equation, 
ie equal to the permeability coefficient K of Carey's equation, sudeti tu- 
tion of 9,678 £ lox for qiL/AD? can be made in equction (9) above ziving 


6 «22678 x 207% 60) 
G1 - 6) g kt a 


Asguming thet the constant k' would be reflected in a recasting of the 
porosity function, @ new expression of this function can be substi tuted 
for ng? /(2 ~ g)* in equation (10) above giving 


(#)" = & / 1.083 ¥ 10°" (22) 
or, by rearrangement, 
Ko wm 1.013 X 10° a (f)” (12) 


where nis a constant applicable only to s particular formation or section 
of a forvation, or & porous nediume 

Table TV, Anpandiz J, has been prepared from caleulations derived 
from the apclication of Fettke's 1934 data to the functions ef equation 
(12) above. Figure 4, Appemiix II, 1e a log-lez plot of these functions, 
x/a®/A.,013 X 10° versus g, which forms a straizht line with a slopa of 4, 


Since 


“f dog t/a°/Aeoh3 x 10° 


then n= &, the slope of the straight line. 
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TV. COMOENT ARO UL SOUSSIGN GF BRSULTS 


Before discussing the results, it should be oointed oat here thet 
the values of porosity data used in this stady represent absolute 
poroeity. However, Fettke” nas shown that in a eajority of cases the 
aifference betwen ahbeclute porosity and a lower value for effective 
percaity in the Draifori Thin sand is but a matter of 0.1 ser cent 
Sverage. There are a fow isolated cases of a deviation alightly greater 
than Goi per sent up to about 1,6 per cent. For purposes of this stady, 
the values of effective perosity sere assumed to be the sane as those for 
absolute poresity, This fact can explain some of the slight deviations 
of the points used to plet the varicus lines an! curves shown by Fiirares 
1, 2, 3 and hk. Other devietions of the points can best be expleined by 
the fact that at the time the data were jenerated, the sccuracy of the 
Séeterninationa of values for permeability, porosity end crein size was not 
nearly as good as that of the present day, due to greatly inproved and 
Standardized sethods for detersining such valuce. 

She regulte of the author's first investiczation clearly show that 
the uverage pero diaaster in a rather [Mineegralzed, consoliteted seilus 
d@ergases consi:iierably from the averase for unconsolidated media, state 
by wutting?? to be reughiy one-fifth of the aversze ;rein diaseter. ‘This 
cen be best exclained by the nen<uniformity of particle size and shape, 
the fine zrain sire, and by the degree of lithi fMesation or censnting 
foand in contolidated welia, It is considered by the author that the 
iatwr exerts the createst effect upon average pore dlenmeter since cerent- 
ing to any extent tenda te restrict, reduce and plug the sore channels of 
the rock. jo utilise this asthe! of correlation of pormenbllity, porosl ty 
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ani crain sise, 1t would be necessary to igtemine the relationship be- 
tween averase effective pore dlmaeter and average grain diameter for wach 
particular reservoir or seotion of the reservoir. 

[n the second ecproseh to the proview, the results of the in- 


veetivation iniicate that the vraine af the Noseny-Carman constant. for 





onsolijated aeiiea ie wueh hizher than the assigned valua of 5, Py way 
of axa sple, vreayuee the consolidated medium to be comoped sencrally of 
sharp particles «here the valae of the ratio of surface to volume factors, 
hg/iyy 18 ToG. 1% nes been previously shown that S$, 2 4,/%,1, Using as 
@ selection fron Fettke's 193i data, d = .093 cantinetar, bh w -0034 arey, 
and @ 2 0145, it is poseible, vy resciution of equation (4) of seation IIT 
of this study, to solve for a value of the Aozeny-Caraan constant in lieu 
of the value of 5 neraally used, %, is resolved as 7/.0093 or 753. 
Seivinsg for the unknown Aoseny-Curnan Gonstent using the above values in 
the equstion gives a value of approxinately 220, by using some of the other 
gets of values for parmaability, porosity, and grein size from the 193, 
data to solve for the hoseny<“araan constant in the save sanner as described 
above jives values for the constant in the zenerel vicinity of 226 wi th 
some variations. The results of similor application cf Fettke's 1933 data 
Were not as conclusive, with the values showing a mich witer variation, 
3%) to 107, for the Lozeny-Caraan constent,. 

43 in the second approsch, the third epproach to the problem also 
imdleated that the Hoseny-Ccarzen constant for consolidated madia ia each 
higher than 5, It hes been previously stated in Section LI! ef this 
atudy that the slope of Line 1 of Figure 3, Appendix II, wae 0s0001052 and 
equivalent to the funetion 1/(2,/ty)"k where k is the Keseny-Carsan 
constant. If, for example, it is asain ansuned that 2,/Z, has a value of 
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7.0 fer this medium, 12 is possible to aclyvs for a value of ether than § 
fer the Koveny~-.araan constant, in this ease, the value wea resolved aes 
194% for the congiant. ince @ small deviation in the Jrawing of line 1 

of either Figure 2 er 3 regalts in «a change ln the value of the slope of 
the line, it can be preeume! that this fact censed the li ffareree in 
values for the /osenyeCarman constant determeined by the tac investigations, 
Therefore, :% can also be presueed thet in the cave of the 17%, data used 
the avere;¢ value for the foxeny-“ersan constant lies sccewhere between 
194 and 220, No siadler esaceese wes achieved in usin: Yettkets 1939 data, 
it 19 believet that the i?3 4ata were net traly representative of the 
reservolr sand nor were they im themselves complete enough to persit érae- 
ing ary conclusions. 

The fourth epproach to the protles was an atteast to show the 
porosity function as a variable rather than 6 constant fsebor in liea of 
a variable for the Koreny~"arsan conetant ea shown by the secon? and third 
investi;sations. An application of the 193; data to equetion (12) of part 
tll dmtlestes thet the mn of the poresity function, #, is epproxinately 
eqaal to 5, as deteratned from the elope of the line slotted in Figure kh, 
Appendix iJ, with higher and lewer variations in some cases, ince the 
perosity fanction is an arbitrary selection by the author, it follows 
that, thie partiealar function, being an exponential, could haw a variety 
of foras all of whieh would, when solved mmerically, cive suvstentialiy 
tha sane values. ¢ waa acestad because of its elaslictty if for ns 
ether reagon. Althouvh this setho!? shews that the perosity function can 
well be a variable factor, it 1s believed thet the worksof Seve am! Bruce”? 
ani ethers have contlusively shown that the porosity function recains 
constant wille the hoseny-Carman constant becomes a variable factor in 
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condolidated media. He attempt wae sade in the fourth lavsstigation to 
make use of fettke's 1938 data because of previous ineonsequentiel 
resul to. 

in ali of the author's four investi cations, it has been show 
that the precedures, equations, and constants applicable ta the empirical 
formulations derived from the studies of unconsolidated porcves media aso 
not heli true in the case of consoalitated perous media. Howewer, it algo 
has been shown that for a liven or particuler censclidseted porous rock 
formation, the eaplrical fcornmulstions apply shen the nornsl constants are 
changed to thome values which fit the siteation, as in the case of thie 
study of the consolidated Uradferd Third send. The evidence toes nat 
conclusively preve whether the porosity function should be a variable er 
a constent factor in peresesbility, nerasi ty and crain size relationship 
formulations, “ueh more inreatiyution ani research cust be carried out 
alon: theaea lines before @ conclusion can be reeehed 2a to which theory 
is correct, although as mentioned in the preesding paragraph present 
eri tence pointe to the poreslty function being a conatant foctor, It is 
very possible thet beth the porosity function ani the Loseny-Carsan ceon- 
giant wicht be found te be verlables is the case of consoliiated porous 
maila., 
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The fleld of snail particle atatistics and maagurevents ia exten- 
sive, complex and relatively unknown. Yor these reasons, the stady of 
particle size im relation to permeability ami porosity of porous media is, 
of iteelf, extensive, complex anid relatively anknown, esceclally in the 
case af conseli:dated porous nedia. fhe suthor avracs with tasient* who 
stated, in effect, that the use of snall pirticle atatiatics ani seasure- 
aents in <etersinations of permeability was so complex an! unresolved that 
che cregently utilized asthois of detersinin: pemmeahility and porosity 
fer cutweich any acthe! using crwin size or ite atteniant preperties to 
ietersine same. However, the author feels thet 14 1a essentially im 
pertent to intensively study the fleld of particle sise in order te zain 
a thorough knowle:e ami avierstanding of macroseepic fMuit flow in porous 
media in the licht of ;rain sise and 1te atteniant factors of packing, 
uniformity of eive and shape, and degres cf lithifiecation, all of which 
affect fluid flew in such eedia, 

Through an intensive study of this nature, 1t is possible that 
some day a ferwalation will be developed whieh will trely hold unter al} 
condi tions and in any circumstances whsrety tho physical proserties of 
permeability, porosity, end crain elise with ite attemiant factors, can be 
euscesefully cerrelated. The four methods of correlation espl«yed in 
thia atudy can be ased in a reetricted degree when applied to the cen- 
solliated porous sedia of particular reservolre or of sections of the 
recarvoira only after considerable reeeareh to Jeterzine the proper 
cangtante apclicable in eech case, In the same rein, and of considerable 
interest, is the method seveloped by tyder”* for aukin; permeshilt oy 
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seasaronenia without the uea of core sanples by means of an eopiricaily 
develope: formula, god only for ome particular sand, which uses only 
Jeternined welvhtga of sieved sani particies which have been passed 
threach a seriece of U. 5, Standard sievas, 

_ Another pesaible investigation or aethod, sugzeated to, but un- 
tried by the author, would be to plot the results of sieving data from 
a core sazmple (per cent by weicht of sanple throegh and on each screen 
sise) versus ‘he average sereen aise for that per cent by weight. Fron 
the resultant carve, selec? a "representative" cer cent «alue and deter~ 
alne the cerreapondin: average screen size or presumed averave grain 
siameter, Using equation (2) of Section III of thie stady, determine 
whieh per eent by weight value an! corresponding screen sise or presumed 
averace crein size can be used to satisfy the equation, ani then detereline 
whether or not. this saee per cent vaiue can be used in eennection with 
all core sazple data from cores taken from the sene formation or section 
of the formation. 

A gumnation of the resalts of this study (oes not reflect the 
original purpove of the study. [4 wes thought thet an evaluation sould 
be meade of the interrelationshipe of the physical prosorties of porous 
mck, permeability, poresity ond grain sise, as well as a sciification of 
existing formulations used te expresa these interrelationahios of phyaies: 
preperties. Ry using date from hishiy ecensolidatet sors wella, 1% was . 
hoped that if a mofifieation sere developed for effectively axpresain; the 
imterrelatilonship of perseability, porosity, ond gruin size, it would hold 
urier any clreumetances ent conmlitions for any type medium, either gore 
solijated or anconsolidatei. However, the comlexities of grain site 
stuty, the iearth of sufficient and eccurate data, and the apparant non 
eonforsity of aprlication of meidifled formiations to consolidated poroun 
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modia data made such endeavors impractical em! Ltmpcasible,. 

Tt is believad that the resulta of this study, although meager, 
can contribute some benefits in a better unlersianding of the problem. 
the ratio of pore volume to grain size diameter hes been shown to decrease 
markedly with an increase in ceventation and a corresponding decrease in 
grein size such as foun? in a consolidated perous media like the Bradford 
third send. ‘The value of the ozenyeCarsan constant hes been shown to 
inereese te very high vélues in consolidated porous media, this is prob- 
ably due to the large increases in values of the tortucstty fector which 
becomes greater as the rock becomes more consolidated because of an in- 
ereage in the effective lenzth of flow path. frobably the fost ianortent 
contribution of this stady ie the realization of the importanee of the 
Koneny~Cargan equction end ite aodifications when apolied te problems of 
fluid flow through porous aciia, Jt holds great possibilities, when sub- 
joct to the proper modi Nication, in ics application to studies of this 
type ani gives promise of eventual sclution of the sany probheme of fluid 
flow still facing the petroleuwn imiustry. 
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irein Sise Permeability Poros ty nf c/o)? Conatant (¢ 2 etn /oy 
d (om) sk (fmretes) ¢ (per cent) xlonm4 (en) x 107% 
(X92 e0GOS 10.2 0678 5 G19 3,82 
oGO52 ° 12,5 -0033 id 3056 
0 CO93 . + 01530 . 3.60 
009% e032 15.3 005190 2278 & 12,80 
20099 eCOL2Z 13k 000395 .09K6 f 32 
20106 eOC3P 15.8 002320 6.1523 " 8,57 
0109 20037 Li 2 O2bb0 1562 ® 78 
eGRL0 08039 1602 eO2hl0 21552 " 8.7% 
oF.20 eOL39 e502 04590 22625 ad Lae75 
eCL21 2G213 2146 009870 3142 ° 17.40 
° 00366 2247 016100 4012 . 22.50 
oOL32 02323 21.7 o1h900 3860 . 21.00 
Oh 35 nO226 Lie? eGL770 1330 . Tels? 
00439 COLT 1560 COLL .1063 . 5.7 
Fettko's 1934 -ata (EK is Caleulatea) 
5619 36k2 
* § 02 
- Fe23 
“ be Th 
" 6.15 
“ 9.20 
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¥ 9-26 
. 15.50 
03225 ° 13.10 
0293 4 23624 ad 
00132 00232, Zhe? 03256 . 15.30 
00135 «0026 1.7 081770 1330 . 7 oly 
0139 0029 15.0 : 339 a 7280 
Fettke's 2933 “ata (K is fxperlmental) 
DOL 000278 1303 002090 =o LS 52519 Jel2 
00042 0007 12,1 000389 «=.052h ’ 3.50 
0 COU2 060131 12,38 eOlhlO 6 i287 ° 6,68 
00042 200550 Lhe3 e030 oI ° 11.00 
o® eiAVL2 15.2 005980 =o 2luhksS sd 13.72 
(Digg 000834 ik.3 eO3740 = 193% 7 10.86 
20087 eO2 752 1507 eATS3O  oht.37 ? 23450 
20051, 15.7 017530 | 4199 ° 23.00 
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Fettke!s 1934 date (K is experimentel) 
a 
drain Sive Permeability Poresity Poa? afm Pas)” , 
d (om) KE (nevreies)  § (ver cant) K 1.013 % lo 
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TABLE LY 
(K is experiaentel) 
Grain Size x /a® jae . Perosity 
d (em) 1.613 x 10 f (per cent) 
0092 5.91 5.83 x 107° 1068 
20092 §.91 §,83 12,5 
«0093 39430 38.70 1b. 
OOD 92. 92.60 15.48 
eCOP9 12, 12.08 130k 
201.06 Bie 3.30 16,8 
eGLO? 31.20 30.75 15.2 
09120 326 31.675 1502 
2ORLO 115.00 113,50 2062 
G21 146,00 148.00 21.6 
0128 223.00 220.00 2247 
00132 185.60 183,00 21.7 
20135 Lise 30 14.10 Lye7 
00139 8,80 8,68 15.0 
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